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Motivated by the experimentally observed y/5 x y/5 iron vacancy order and a block spin antiferromagnetic 
phase with large magnetic moment in K0.sFe1.6Se2, we study the magnetic phase diagram of an extended 
J 1 — J 2 model on a | -depleted square lattice with y/E x y/E vacancy order, using a classical Monte Carlo 
analysis. The magnetic phase diagram involves various antiferromagnetically ordered phases, and most of 
them have higher order commensuration. We find that the experimentally relevant block-spin state occupies a 
significant portion of the phase diagram, and we discuss the spin dynamics of this phase using a linear spin- 
wave analysis. By comparing the calculated magnetization with the experimental values of magnetic moment, 
we determine the physical parameter regimes corresponding to the block spin antiferromagnetic phase. Based 
on our spin wave calculations in different parameter regimes, we show how spin-wave degeneracy along the high 
symmetry directions of the magnetic Brillouin zone can provide information regarding the underlying exchange 
couplings. We have also analyzed the magnetic phase diagram of a Ji — J 2 model on two different modulated 
square lattices relevant for K !/ Fei.5Se2, which respectively exhibit ^-depleted 2x2 and 4x2 vacancy ordering. 



I. INTRODUCTION 

Ever since the discovery of the high temperature supercon- 
ductivity in iron pnictides in the vicinity of an antiferromag- 
netically ordered phased, the strength of electronic correla- 
tions in iron based superconductors has been a central issue 
of debate. Significant efforts based on both weak and strong 
coupling pictures have been made to explain the (tt, 0) mag- 
netic order in the parent compounds, and the emergence of su- 
perconductivity upon carrier doping. On one hand, the (tt, 0) 
collinear antiferromagnetic order arises in a weak coupling 
picture that relies on the presence of quasi-nested electron and 
hole pockets respectively at M and T points of the extended 
Brillouin zone—. On the other hand, the experimentally ob- 
served large electrical resistivity ("bad metal" with kpl ~ 1), 
a strong suppression of Drude weight 4 ., and the temperature- 
induced spectral weight transfer— ~— have suggested sizable 
strength of electronic correlations. These experimental evi- 
dences have suggested the placement of the iron pnictides in 
close proximity of a putative Mott transition 8 - - — . In a metal- 
lic system close to a Mott transition, quasi-local moments are 
expected to arise, which for the iron pnictides are described 
in terms of J\ — J2 couplings^; here J\ and J2 refer to 
nearest-neighbor and next-nearest-neighbor spin exchange in- 
teractions on the iron square lattice. Hence in the strong cou- 
pling scenario, the (71", 0) collinear antiferromagnetic order in 
these materials can be described in terms of such a J% — J2 
model. This picture is further supported by the experimen- 
tal observation of zone boundary spin wave excitations in the 
magnetically ordered state at low temperatures 11 . When pro- 
tected by the broken symmetry, many low energy properties 
of the ordered state have been addressed using either strong 
or weak coupling approaches, and this has contributed to the 
dilemma over the strength of electronic correlations. 

For 11-iron chalcogenides FeTei-^Se^, a nesting based 
approach fails to describe the magnetic ground state— — . 
These materials have been shown to exhibit both commensu- 
rate and incommensurate antiferromagnetic state with a large 
magnetic moment ~ 2.0^b- However it should be noted that 



in spite of the stronger correlations in 1 1 -chalcogenides in 
comparison to pnictides, the parent compound FeTe is still 
metallic. 

The discovery of superconductivity in the new 122 iron- 
chalcogenides (K, Tl, Cs, Rb^FezSe^— , in the vicinity 
of an antiferromagnetic state with aT c ~ 30A', promises to 
shed new light on both the role of electronic correlations and 
nature of magnetism and superconductivity. In these 122 iron- 
chalcogenides, there are parent compounds which are antifer- 
romagnetically ordered and insulatin g 18 21 . Furthermore, both 
the angle resolved photoemission experiments 2 ^ 4 - and band 
structure calculations^— reveal unique fermiology. Unlike 
the pnictides, there are no hole Fermi pockets; the Fermi sur- 
face is entirely made up of electron Fermi pockets, with the 
dominant pieces located near the M points of the extended 
Brillouin zone. The bare magnetic susceptibility xo(q, w) for 
these new materials, as determined from the band structure, 
is weakly peaked around q = (tt, tt), and can not explain the 
origin of the observed complex magnetic order with a large 
magnetic moment ranging from 2/j.b to ~ 3.4/is which is 
intertwined with iron vacancy order— — . The observation of 
an insulating state and the emergence of antiferromagnetic or- 
der with large moment in the absence of quasi-nested electron 
and hole pockets certainly point towards stronger electronic 
correlation effects. The possible band narrowing effects due 
to vacancy ordering, and a resulting Mott insulating phase for 
the parent materials has been discusse d 32,33 . 

In this paper, we study the magnetic order and dynamics of 
an extended 3\ — J2 model on modulated square lattices. For a 
lattice (called L\ below) appropriate 18 for the vacancy order- 
ing of a K l/ Fei.gSe2, our analysis builds on the important in- 
sights introduced by Cao and Dai 34 that vacancy ordering may 
modulate the exchange interactions. We show that the experi- 
mentally observed block-spin stats 2 ^ 3 -^ occupies a significant 
portion of the phase diagram, and we discuss the spin dynam- 
ics of this phase using a linear spin-wave analysis. The large 
unit cell with multiple iron ions leads to both acoustic spin 
wave and gapped optical spin wave modes. Based on our spin 
wave calculations in different parameter regimes correspond- 
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ing to the block spin antiferromagnetic phase, we show how 
the measurements of spin gaps and degeneracy along the high 
symmetry directions of the magnetic Brillouin zone provide 
valuable information regarding the signs of the underlying ex- 
change couplings. We have also analyzed the magnetic phase 
diagram of a J\ — J2 model on two different modulated square 
lattices (called L2 and L3 below) potentially relevan t 18 ! 35 to 
Kj / Fei.5Se2, which respectively exhibit ^--depleted 2x2 and 
4x2 vacancy ordering. We note that we have not considered 
the further-neighbor coupling J 3 , since the J\ — J 2 couplings 
appear to be adequate for K a Fe x Se2. 

Our paper is organized as follows. In Sec. II we introduce 
an extended J\ — J2 model, as appropriate for the modulated 
square lattice L\ with Vo x y5 vacancy order (Fig. |l(a)| i. 
In Subsec. II A and Subsec. II B we respectively discuss 
possible magnetic phases and phase diagram of the extended 
J\ — J2 model for the L\ lattice, based on a Monte Carlo 
calculation using classical spins. In Sec. Ill we discuss the 
spin wave results for the block spin antiferromagnetic phase. 
In Sec. IV we introduce the appropriate J\ — J2 models for 
the modulated lattices L2 (Fig. |6(a)| i and L3 (Fig. |6(b)[ i. We 
discuss possible magnetic phases and phase diagrams for the 
L2 and L3 lattices, respectively in Subsec. IV A and Subsec. 
IV B. A summary is given in Sec. V 

H. EXTENDED Ji - J 2 MODEL FOR Li LATTICE WITH 
\/5xV^ VACANCY ORDER 

The modulated lattice L\ with y/E x y/5 vacancy order is 
shown in Fig. 1 1 (a)| Shown here is the structural unit cell com- 
prising four Fc sites and a larger unit cell consisting of eight 
Fe sites. The eight site unit cell is the magnetic unit cell of 
the experimentally relevant block-spin antiferromagnetic state 
and this also turns out to be convenient to describe other an- 
itiferromagnetic phases to be specified below. In Fig. |l(b)| 
we show the structural Brillouin zone (SBZ) corresponding 
to the four-site unit cell and the magnetic Brillouin zone of 
the block-spin antiferromagnetic state (MBZ1) correspond- 
ing to the eight site unit cell. A recent DFT calculation 36 
shows that the Fe vacancy order induces a tetramer lattice dis- 
tortion which may strongly influence the superexchange cou- 
plings. This distortion makes the intra-block and inter-block 
exchange couplings different, and leads to the following mod- 
ulated Ji — J2 Hamiltonian^i 

Hi = — I Jl E ■ Sip + J 2 ^2 Sia ■ Si/3 i 

i 1 (ccfS) «a/3)> J 

+ \ E { J i E s <° • s ^ + J2 E s <° • s i4' (1) 

ijtj L (a/3) {(a/3)) J 

where the Latin indices i, j correspond to the unit cells, each 
enclosing a single vacancy and comprising four Fc sites. The 
Greek indices a,/3 = 1, 2, 3, 4 correspond to the four Fc sites 
of the unit cell. We have respectively denoted the nearest 
and next nearest neighbors by (a/3) and ((a/3)). Within the 
same unit cell, the exchange couplings between the nearest 
and next nearest neighbor iron sites are respectively denoted 



by Ji and J 2- Between the neighboring unit cells, the cou- 
plings between the nearest and next nearest neighbor iron sites 
are respectively J{ and J' 2 . In order to incorporate the modula- 
tions of the exchange couplings by the vacancy order, we will 
allow the exchange couplings to take both antiferromagnetic 
( J a , J' a > 0) and ferromagnetic (J a , J' a < 0) signs. 



A. Magnetic phases for L\ lattice 

To determine the phase diagram of the extended J1-J2 
model defined by Eq.Q] we apply classical Monte Carlo tech- 
nique. We treat the spins as classical 0(3) vectors, and adapt 
the standard Metropolis algorithm. The Monte Carlo simula- 
tions have been performed on finite lattices with linear dimen- 
sion up to 80, and at temperatures as low as T ~ 10~ 3 J' 2 . By 
checking both the static spin structure factor and the Monte 
Carlo snapshots of spins in real space, we identify the order- 
ing wave vectors and corresponding real-space spin pattern of 
various phases stabilized at low temperatures. We have found 
seven possible magnetic phases. One of these is an incom- 
mensurate magnetic phase (no higher order commensuration), 
with very broad structure factor in the entire EBZ, SBZ and 
MBZ1. This phase occupies a small sliver of the parameter 
space in the phase diagram and will not be discussed in detail. 
In the magnetic phase diagram we denote this incommensu- 
rate phase as IC. The other six magnetically ordered phases 
have commensurate ordering vectors in the MBZ1. These 
phases are denoted as AF1, AF2, AF3, AF4, AF5 and FM in 
the following discussions, and are shown in Fig. [2(a)] through 
Fig. |2(f)| The AF1 phase corresponds to the experimentally 
relevant block-spin antiferromagnetic state. 

All the six magnetic states shown in Fig. 2 can be described 
by the following single-mode spiral: Si a = a;cos(Q • R^ + 
@a) + V sin(Q Hi + 9 a ), where Rj and Q are respectively the 
position vector in the 8-site unit cell, and the ordering wave 
vector in the MBZ1 . The canting angle at site a is denoted by 
9 a . The fact that these six states are the lowest energy states, 
has also been confirmed by minimizing the ground-state en- 
ergy with respect to Q and 9 a . The commensurate ordering 
wave vectors, the canting angles and the classical ground state 
energy per site are listed in TABLE. I. 



The four phases AF1, AF3, AF4 and FM are described 
by the same ordering wavevector Q = (0,0) in the MBZ1, 
whereas Q = (71", 7r) for AF2 and AF5. The phases with 
same ordering vectors are distinguished by the assignment 
of the canting angle values. Depending on the canting an- 
gle values, we obtain collinear and non-collinear magnetic 
states. The AF1 phase is the collinear, block-spin antiferro- 
magnetic phase found in the neutron diffraction experiments 
on Ko.8Fei.6Se 2 . On a four site unit cell all the spins are fer- 
romagnetically aligned and form a large block spin, and these 
block spins on the adjacent four site unit cells are antiferro- 
magnetically aligned. This state can be equivalently described 
as a (71", 7r) AF ordered phase in the SBZ. In the following sec- 
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FIG. 1. In panel (a) we show the L\ lattice with y/E x \/5 vacancy 
order, with the shaded circles describing the vacancies. The struc- 
tural unit cell due to the vacancy ordering is marked by dotted green 
lines, with four Fe sites residing on a plaquette 1234 which is marked 
by thick, dark lines. The intra unit cell exchange couplings Ji, J2 
and inter unit cell exchange coupling J[ , J2 are also illustrated. The 
larger unit cell consisting of eight sites 1, 2, 8, and marked by 
dashed green lines corresponds to the magnetic unit cell of the AFl 
phase. In panel (b) we show three different Brillouin zones (BZ). 
The extended BZ (EBZ) — n < k x < n, —n < k x < tt, marked 
by solid black lines corresponds to 122-unit cell that encloses one Fe 
site. The structural BZ (SBZ) of the modulated lattice corresponding 
to the unit cell with four sites is marked by the dashed blue lines. 
This BZ encloses |-th of the area enclosed by the EBZ. The smallest 
BZ also marked by solid black lines corresponds to the magnetic BZ 
of the AFl phase (MBZ1). The high-symmetry points in the SBZ 
and MBZ1 have the following coordinates: in SBZ, X = |), 
Y = (-f , 2*), M = (f , f ); and in MBZ1, X' = 
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FIG. 2. In panels (a) through (f) we show the real space arrange- 
ments of the spins for six commensurate magnetic states. AFl, 
AF3, AF4 are collinear antiferromagnetic states with ordering vec- 
tor Q = (0,0) in the MBZ1, and these phases are distinguished 
by the different canting angle assignments. AF2 and AF5 are non- 
collinear antiferromagnetic states, with ordering vector Q = (tt, it) 
in the MBZ1. The conventional ferromagnetic state is denoted as FM 
and has Q = (0,0). AFl phase corresponds to the experimentally 
relevant block spin antiferromagnetic phase. 



tion we address the role of quantum fluctuations in this phase 
using linear spin wave analysis. 

The AF2 phase is a non-collinear antiferromagnetic state. 
On each four site unit cell, nearest neighbor spins are or- 
thogonal to each other, whereas the next nearest neighbor 
spins are antiparallel to each other. The canting angle for all 
the spins on the adjacent unit cells are shifted by an angle 
arctan(2 Jj/ J'\)- In terms of the SBZ, the ordering wavevec- 
tor is at both (tt, 0) and (0, tt). Note that this state is different 
from the collinear (tt, 0)/(0, tt) state in the J\—J<x model when 
J2 > Ji/2. The AF2 state does not break the C4 rotational 
symmetry, and it cannot be described by a single ordering 
wavevector in the SBZ. The AF3 phase is a collinear magnetic 
state. On a four site unit cell the spins are antiferromagneti- 
cally aligned and each unit cell has the same spin arrangement. 
Therefore in the SBZ this state has ordering vector (0,0). In 
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Phase 


Q-vector 


9 a 


Energy per site 


AF1 


Q = (o,o) 


$1,2,3,4 = 
$5,6,7,8 = T 


Ji - J'2 + ^ 


AF2 


Q = (tt.tt) 


a (i— l)ir 

_ (i— S)vr 

tan- (f ) 




AF3 


Q = (0,0) 


$1,3,5,7 = 
$2,4,6,8 = TT 


J2+ 2 J[ Ji J* 


AF4 


Q = (0,0) 


$1,3,6,8 = 
$2,4,5,7 = 7T 


^ - Ji + j 2 


AF5 


Q = (7T,7r) 


$1,3 = 0, $6,8 = ■§ 
$2,4 = TV + $0 

$ 5 , 7 = rr/2 + $ () 

0° = tan_1 (4) 




FM 


Q = (0,0) 


6i = for 
i = l,...,8 


Ji + + ^ 



TABLE I. The ordering wave vector Q in MBZ1, the canting angles 
and the ground-state energies per Fe site of the six commensurate 
magnetic states shown in Fig, |2(a)| through Fig, |2(f)| 



the AF4 phase every neighboring spin is antiferromagnetically 
aligned and this phase corresponds to conventional Neel state 
on the 122 unit cell. Also notice that on each four site unit 
cell spins are antiferromagnetically aligned, but on the adja- 
cent unit cells all the spins are flipped. Therefore in terms of 
SBZ AF4 has (it, it) ordering vector. The AF5 phase exhibits 
non-collinear magnetic order. On each four site unit cell the 
canting between nearest neighbor spins depends on the ratio 
J[/2Ji, whereas the next nearest neighbors are ferromagnet- 
ically aligned. Going to an adjacent unit cell all the canting 
angles are shifted by tt/2. Finally the FM phase corresponds 
to the conventional ferromagnetic state. In the following sub- 
section we describe the associated magnetic phase diagram. 

B. Phase diagram for Li lattice 

The seven phases mentioned above give rise to a complex 
phase diagram, and most of the phase diagram excluding the 
incommensurate IC phase, can be obtained by comparing the 
tabulated energies for the six magnetic phases. We measure 
energies in the units of J 2 and in Fig. |3(a)| through Fig. |3(d)| 
we plot the phase diagrams in the J\j ' J 2 — J2/J2 plane, for 
four different values of J[/J 2 ratio. The phase diagram in 
Fig. [3(a)] corresponds to J[/ J 2 = 1 and qualitatively remains 
applicable for J{/ J 2 < 2. The phase diagram in Fig. |3(b)| 
corresponds to J[/J 2 = 4 and is qualitatively applicable 
for any J[/J 2 > 2. The phase diagram in Fig. [3(c)] cor- 
responds to J'l/ J 2 = —1 and is qualitatively applicable for 
any J[/J 2 > —2. The phase diagram in Fig. [3(d)] corre- 
sponds to J[/J 2 = —4 and is qualitatively applicable for any 
J'l I J'2 < — 2. Considering the symmetry of the Hamiltonian 
in Eq. [TJ we find the ground states for J[ < or J 2 < can 
be obtained from the states in the phase diagram at J[ > 
and J 2 > by performing the following two transformations: 



a) sending S Q — s- — S Q for a = 2,4,5,7 in the 8-site unit 
cell, and J\ — > —J\, J[ — > —J[\ b) sending S a — > —S a for 
a = 2,4,6,8, and J\ — > —J\, J 2 — > —J 2 - The experimen- 
tally relevant block spin state AF1 occupies a large portion of 
the phase diagram. Only for J[/J 2 < —2 is the AF1 phase 
absent. This provides a simple constraint on the J[/J 2 ratio. 
It is important to note that AF1 phase extends into a small re- 
gion where all the exchange couplings are antiferromagnetic. 
We further note that the values of the exchange couplings pre- 
dicted in Ref. [34] based on LDA calculations are close to the 
AF1-AF2 phase boundary. This suggests, on one hand, it is 
possible that in certain materials J 2 can be tuned to cross the 
AF1-AF2 phase boundary to stabilize the AF2 ground state, 
which can still be described by our theory. On the other hand, 
it also suggests considerable quantum fluctuations in the ex- 
perimentally observed block-spin (AF1) state. In the next sec- 
tion we investigate the quantum fluctuations in the AF1 state 
within the framework of linear spin-wave theory. We demon- 
strate how the spin wave dispersions in the AF1 phase depend 
on the underlying exchange constants in the three regimes cor- 
responding to Fig. |3(a)| through Fig. |3(c)| 




1/1; J / J ' 2 
(c) (d) 

FIG. 3. In panels (a) through (d) we respectively show the mag- 
netic phase diagrams in the J1/J2 — J2/J2 plane for J1/J2 = 
1, 4, —1, —4. The phase diagrams in panels (a), (b), (c) and (d) are 
respectively applicable for J1/J2 < 2, J[/ ' J 2 > 2, J[/ 1 J' 2 > -2, 
J'l I J'2 < ~ 2, with only quantitative shifts of the phase boundaries. 



III. SPIN WAVE RESULTS FOR BLOCK SPIN 
ANTIFERROMAGNETIC STATE AF1 

To study the effects of quantum fluctuations on the block- 
spin antiferromagnetic state using linear spin-wave theory, 
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we introduce the following linearized Holstein-Primakoff 
(HP) transformations. Within the 8-site unit cell, we have 

Sf a = S - a] a a ta , 5+ = \/2Sa ia , S~ a = V^Sa\ a 
for a = 1,2,3,4, and Sf a = -S + a\ a a la , S+ = 
V2Sa\ a , S~ a = V2Sa la for a = 5,6, 7, 8. After per- 
forming the Fourier transformation, and defining the spinor 

= ( a l,k> a 2,ki a 3,k' a 4,k> a 5,-k!a6,-k,a7,-k, Cl8,-k), the 

Hamiltonian H = ^ k V'k^kV'k can be diagonalized via a 
Bogoliubov transformation 6^, k = X)s=i 2,3 4 U a f3M a \^ + 
V r Q /j,ka4+/3,-k, where ]C/3(I^W,k| 2 - |E4/3,k| 2 ) = 1. 




-2-101 -2-101 

J,/J\ 

(a) (b) 

FIG. 4. Contour maps of magnetization m (in units of /is) in spin- 
wave calculations for J{ = J 2 (in (a)) and J{ = — J 2 (in (b)). The 
dashed green lines are the contour lines of m = 2.0/is per Fe (upper 
line) and m — 3.4/is per Fe (lower line). The region in between 
defines the physical parameter regime. 

It is worth discussing which regime in the parameter space 
is most relevant to the experimentally observed block-spin 
state. To address this issue we may define a physical pa- 
rameter regime in the AF1 phase by requiring the theoret- 
ically calculated magnetization, m, to be within the range 
of magnetic moments determined in experiments. Experi- 
mentally, it is found that the magnetic moments take val- 
ues between 2{1b and 3.4/j,b for various A y Fei.eSe2 {A = 
K,T1, Cs, Rb) compounds. 31 Theoretically, m = S — 5m, 
where 8m = § E ? ^=i,2,3 4 JkeMBZ l^9,k| 2 , is the correc- 
tion to the magnetization due to quantum fluctuations. We 
show the contour maps of m from spin-wave calculations for 
J[ = J ' 2 and J[ = -J 2 in Fig. |4(a)| and Fig. [4(b)] respec- 
tively. In both cases, we find the physical parameter regime is 
only limited within a narrow region near the boundary of the 
AF1 state. Note that the LDA calculation also suggests the ex- 
change coupling parameters are close to the phase boundary 
between AF1 and AF2 states^ 

To discuss the spin-wave spectra, we further follow the 
suggestion by LDA calculation^ to assume both J2 and J 2 
to be antiferromagnetic, and show in Fig. |5(a)| to Fig. |5(c)| 
the calculated spin-wave dispersions for three representative 
points in the phase diagram. The spin-wave dispersions are 
obtained from the eigenenergies in MBZ1. Due to the four 
Fe sites within a block, in addition to the gapless acoustic 
branch (Goldstone mode), there are also three gapped opti- 
cal branches. Each of these four branches are also doubly 




X' y r X' M" r 



(c) 

FIG. 5. Spin- wave spectra along the high-symmetry directions of the 
MBZ1 of the block spin state for three representative sets of param- 
eters: a) Ji = O.IJ^, J 2 = 0.05J^, J[ = J' 2 ; b) J 1 = -0.5 J^, 
J 2 = 0.5J 2 , J{ = J 2 ; c) Ji = -1.5J 2 , J 2 = 0.2J 2 , J[ = -J 2 . In 
all three cases, J 2 = 1, is the energy unit. 



degenerate. Working out the eigenenergies at V and M' points 
in MBZ1, we find that the degeneracy of the eigenenergies at 
these two points is helpful to determining the sign and relative 
strength of the exchange couplings. As shown in Fig. |5(a)| 
when Ji and J[ are both positive, the top two optical branches 
at r point are degenerate, and the lowest optical branch and 
the acoustic branch at M' point are degenerate. In this case, 
there is a finite gap between the top two optical branches and 
the rest of the spin- wave spectrum. If J\ < and J2 ~ |Ji|, 
as shown in Fig. |5(b)| the lowest two optical branches at T 
point are degenerate, and the lowest optical branch and the 
acoustic branch at M' point are also degenerate. In this case, 
a finite gap separates the top optical branches from the rest of 
the spin-wave dispersion. On the other hand, when J\ < 
but J2 <C I Ji|, the lowest two optical branches at both V and 
M' points are degenerate. In this case, the acoustic branch is 
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completely separated from the optical ones. Moreover, a finite 
gap also separates the top optical branches from the two lower 
optical branches, as shown in Fig. |5(c)| 

The spin-wave spectrum for Rb y Fe! 6 Se 2 has been recently 
measured through neutron scattering experiment^ In the ex- 
perimental spin-wave dispersion the acoustic branch is well 
separated from the optical ones, is similar to the one shown in 
Fig. |5(c)| The dispersion is well fitted to the extended J1-J2 
model in Eq.Q]with J\ < and J2 <C |JiJ. This confirms the 
direct relevance of our results to the experiments. 



IV. Ji J 2 MODEL FOR 2x2 VACANCY ORDERED L 2 
LATTICE AND 4x2 VACANCY ORDERED L 3 LATTICE 
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FIG. 6. In panels (a) and (b) we respectively show 2x2 vacancy 
ordered L2 lattice and 4x2 vacancy ordered L3 lattice. In panel 
(a) the magnetic unit cell for L2 lattice is marked by dashed green 
line, and the magnetic unit cell consists of three Fe sites denoted as 
1, 2, 3. In panel (b) the three-Fe magnetic unit cell for L3 lattice is 
shown using dashed green lines, and we also show a larger unit cell 
that consists of six Fe sites denoted as 1, 2, 6. 



In this section we discuss the magnetic phase diagrams for 
two different modulated lattices L2 and L3 which respectively 
possess 2x2 and 4x2 vacancy orders. These two vacancy 
orders have been claimed to be important for K 2/ Fei.5Se2. In 
Fig. |6(a)| and Fig. |6(b)| we respectively show L2 and L3 lat- 
tices. We note that the 2x2 vacancy ordering preserves the 
Ca v symmetry of the square lattice, and for this reason we 
analyze an isotropic J\ — J2 model for L2 lattice. However, 
4x2 vacancy order breaks the C± v symmetry, and for this 
reason we analyze an anisotropic J\ x — J\ y — J2 model for 
L3 lattice. The explicit form of the Hamiltonian is given by 

H2 = Jix ]^ Si • Sj + x + Ji y Sj • Si + y 

i i 

+ ^X) S *' S i+x±V' (2) 

i 

where J\ x = J\ y for L2. For the description of different mag- 
netic states on L2 we use a 2 x 2 unit cell including three Fe 
sites, as shown in Fig. |6(a)| with dashed green lines. To de- 
scribe the magnetic states on L3, we can use a similar three-Fe 
unit cell shown as shown in Fig. |6(b)| with dashed green lines. 
Sometimes it is more convenient to understand the nature of 



the magnetic state on L3 by using a larger 4x2 six-Fe unit 
cell as shown in Fig. |6(b)| with dotted green line. 



A. Magnetic phases and phase diagram for L2 lattice 

We again employ a classical Monte Carlo technique to un- 
derstand the phase diagram for L2 lattice. There are three 
phases Fl, F2, FM and these are respectively shown in 
Fig. |7(a)| Fig. |7(b)| and Fig. |7(c)| The phase diagram as a 
function of the ratio J\ / J2 is shown in Fig. |7(d)| These three 
phases are commensurate in terms of three-Fe magnetic unit 
cell, and all three phases have ordering vector Q = (0, 0). All 
three states can be described by a single mode spiral ansatz, as 
introduced for L\. The values of Q, the canting angles 9 a , and 
energy per unit site for all three phases are listed in TABLE. [XT] 
For general J1/J2 ratio, the Fl state is a non-collinear ferri- 
magnetic state. But it turns out to be a non-collinear antifer- 
romagnetic state if Ji = 2 J 2- F2 state has the conventional 
Neel arrangement, but the presence of vacancy gives rise to 
nonzero magnetic moment per unit cell. For this reason, F2 
is a collinear ferrimagnetic phase. Finally FM is the conven- 
tional, collinear ferromagnetic state. The phase boundaries 
are found by using the energy values listed TABLE. [II] 



Fl 



\ r A ■ r \ 



V 



V 



V 



V 



F2 



X \AA \A 



(a) 



V A A A A 

(b) 



FM 



\ A A AAA 



FM Fl 



F2 



J/J, 



(c) 



(d) 



FIG. 7. In panels (a), (b), (c) we respectively show three possible 
phases Fl, F2, F3 found for the L2 lattice. Fl and F2 are respec- 
tively non-collinear and collinear ferrimagnetic states, and F3 is the 
conventional ferromagnetic state. In panel (d) we show the relevant 
phase diagram as a function of the ratio J1/J2. 
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Phase 


Q-vector 


9 a 


Energy per site 


FM 


Q = (0,0) 


61,2.3 = 


4(Ji + J 2 )/3 


Fl 


Q = (0,0) 


9j = cos-^-Ji^Ja) 

02 =-01,03=0 


-(8J| + J?)/6 


F2 


Q = (0,0) 


01,2 = 7T, 3 = 


4(J 2 - Ji)/3 



TABLE II. The ordering vector Q in the MBZ for three Fe site unit 
cell, the canting angles and the ground-state energies per Fe site for 
the magnetic states shown in Fig.|7(a)|through Fig.|7(c)| 



B. Magnetic phases and phase diagram for L3 lattice 

The possible magnetic states for the L3 lattice are again 
found from a classical Monte Carlo calculation. There are six 
phases denoted as I through VI, and the spin arrangements 
in these phases are shown in Fig. |8(a)| through Fig. |8(h)[ In 
the MBZ corresponding to three-Fe unit cell the spin states I 
through VI have either (0,0) or (0, 7r) wavevector. Interest- 
ingly, we find in phases V and VI the pairs of states (V,V) 
and (VI, VI') are energetically degenerate, and coexist in the 
phase diagram. In Fig. |8(i)| we show the phase diagram in the 
J\xl J2 — J\yl ' J2 plane. The wavevectors, canting angles and 
energy per site for states I to VI are displayed in TABLE. [HI] 

The magnetic states I, II, III and IV have collinear spin ar- 
rangements, whereas V, V',VI, VI' have non-collinear spin 
arrangements. The state I is an antiferromagnet with the spins 
having conventional Neel arrangement. The state II is the con- 
ventional ferromagnetic state. The state III has conventional 
(it, 0) antiferromagnetic arrangement, and even in the pres- 
ence of the vacancy order, this state remains antiferromag- 
netic. To see this one needs to consider the large 4x2 unit 
cell with six Fe sites. The state IV has conventional (0, ir) an- 
tiferromagnetic arrangement. But, the vacancy order leads to 
nonzero magnetic moment for 4x2 unit cell, and this state 
becomes ferrimagnetic. Interestingly, state IV and state III 
are energetically degenerate when J\ x = J\ y , though the L'i 
lattice breaks the C4 symmetry. This degeneracy is not lifted 
even by introducing a finite J3 coupling between the 3rd near- 
est neighbor spins. The state V is a non-collinear ferrimag- 
netic state, and state VI is a non-collinear antiferromagnetic 
state. Both of them can be described by the single mode spiral 
ansatz and have commensurate ordering wavevectors in the 
three-Fe MBZ. But states V and VI' cannot be described by 
the single mode spiral ansatz in the three-Fe MBZ, though 
they are degenerate with V and VI, respectively. By compar- 
ing Fig. [8(e)] with Fig. 8(g) we see that state V can be obtained 
from V by shifting all spins in the An + 1 rows to the left by 
one lattice spacing. Similarly, state VI' can be obtained from 
VI by first flipping all spins in the 4n + 1 rows and then shift- 
ing them to the left by one lattice spacing. Note that both state 

V and state VI' can be described by the single mode spiral 
ansatz and have commensurate ordering wavevector (0, 0) in 
the BZ corresponding to the 4x2 six-Fe unit cell. Both state 

V and state VI' are non-collinear ferrimagnetic. 
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FIG. 8. In panels (a) through (h) we show the spin arrangements for 
possible magnetic states I through VT realized for L3 lattice. States 
I, II, IB, and IV are collinear and the states V, V, VI, VP are non- 
collinear. The states IV and V are ferrimagnetic. The states I, III, 
and VI are antiferromagnetic, and the state II is ferromagnetic. The 
pairs of states (V, V) and (VI, VP) are degenerate, and coexist in the 
phase diagram displayed in panel (i). 
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Phase 


Q-vector 


9 a 


Energy per site 


I 


Q = (o,jt) 


81 = tt, 6 2 ,z = 


2(2 J2 — Ji x — Jly)/3 


II 


Q = (0,0) 


#1,2,3 = 


2(J lx + J ly + 2J 2 )/3 


TTT 
111 


O (C\ *7r\ 

14 — (V, TV) 


"2 — 7T 5 "l,3 — U 




IV 


Q = (0,0) 


6*1.2 = TV, 3 = 


2{J\ X — Ji y — 2J 2 )/3 


V 


Q= (0,0) 


#2 = #3=0 


-[8Jl + (2J 2 + 


VI 


Q = (0,jt) 


#2 = 7T - 01, #3 = 


[8JL + (2J 2 - 

Jly) ]/l2Ji x 



TABLE III. The ordering wavevector Q in three-Fe MBZ, the canting 
angles, and the ground-state energies per Fe site for the magnetic 
states shown in Fig.|8(a)|through Fig.|8(h)| 



V. SUMMARY AND CONCLUSION 

To summarize, we have studied the magnetic phase diagram 
of an extended J\ — J2 model on several modulated square 
lattices. For a ^-depleted square lattice with y/E x y/5 va- 
cancy order we have shown that a block-spin antiferromag- 
netic state, of the type observed in Kj,Fei.6Se2, arises over 



a significant region of the phase diagram. This region com- 
prises three parts: a) all the intra-block exchange couplings 
Ui,J2) and inter-block couplings (JJ,./ 7 ,) are antiferromag- 
netic; b) J 2, J[,J^ are antiferromagnetic, while J\ is ferro- 
magnetic; c) Jj and J 2 are antiferromagnetic, while Ji and 
J 2 are ferromagnetic. 

We have also calculated the spin-wave spectrum and the 
renormalized magnetic moment in the block-spin state, from 
which we locate the most experimentally relevant parame- 
ter regime in the phase diagram. By studying the spin-wave 
spectrum in this experimentally relevant parameter regime, it 
is suggested that measurements of spin gaps and degeneracy 
along the high symmetry directions of the magnetic Brillouin 
zone will provide valuable information regarding which part 
of the parameter space the exchange couplings belong to. 

Finally, we have studied the magnetic phase diagram of a 
J\ — J2 model on j -depleted square lattices with 2 x 2 or 
4x2 vacancy orders. These phase diagrams are likely to be 
relevant for K y Fei.5Se2. 

This work was supported by NSF Grant No. DMR- 
1006985 and the Robert A. Welch Foundation Grant No. C- 
1411. During the final stage of writing up this paper, related 
works on the magnetism of K a Fei.6Se2 have appeared^^S. 



1 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. 
Chem, Soc, 130, 3296 (2008). 

2 Z.-A. Ren, W. Lu, J. Yang, W. Yi, X.-L. Shen, C. Zheng, G.-C. 
Che, X.-L. Dong, L.-L. Sun, F. Zhou, and Z.-X. Zhao, Chin. Phys. 
Lett., 25, 2215 (2008). 

3 C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. Ratcliff II, J. L. 
Zarestky, H. A. Mook, G. F. Chen, J. L. Luo, N. L. Wang, and P. 
Dai, Nature 453, 899 (2008). 

4 M. M. Qazilbash, J. J. Hamlin, R. E. Baumbach, L. Zhang, D. J. 
Singh, M. B. Maple, and D. N. Basov, Nat. Phys. 5, 647 (2009). 

5 W. Z. Hu, J. Dong, G. Li, Z. Li, P. Zheng, G. F. Chen, J. L. Luo, 
and N. L. Wang , Phys. Rev. Lett. 101, 257005 (2008). 

6 A. V. Boris, N. N. Kovaleva, S. S. A. Seo, J. S. Kim, P. Popovich, 
Y. Matiks, R. K. Kremer, and B. Keimer, Phys. Rev. Lett. 102, 
027001 (2009). 

7 J. Yang, D. Huvonen, U. Nagel, T. Room, N. Ni, P. C. Canfield, 
S. L. Bud'ko, J. P. Carbotte, and T. Timusk, Phys. Rev. Lett. 102, 
187003 (2009). 

8 Q. Si, and E. Abrahams, Phys. Rev. Lett. 101, 076401 (2008). 

9 J. Dai, Q. Si, J.-X. Zhu, and E. Abrahams, Proc. Natl. Acad. Sci. 
106, 4118 (2009); Q. Si, E. Abrahams, J. Dai, and J.-X. Zhu, New 
J. Phys. 11, 045001 (2009). 

10 A. Kutepov, K. Haule, S. Y. Savrasov, and G. Kotliar, Phys. Rev. 
B 82, 045105 (2010). 

11 J. Zhao, D. T. Adroja, D.-X. Yao, R. Bewley, S. Li, X. F. Wang, 
G. Wu, X. H. Chen, J. Hu, and P. Dai, Nat. Phys., 5, 555 (2009). 

12 J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. Wu, G. F. 
Chen, X. Dai, J. L. Luo, Z. Fang, N. L. Wang, EuroPhys. Lett. 83, 
27006 (2008). 

13 F.-C. Hsu, J.-Y. Luo, K.-W. Yeh, T.-K. Chen, T.-W. Huang, P. M. 
Wu, Y.-C. Lee, Y.-L. Huang, Y.-Y. Chu, D.-C. Yan and, M.-K. 
Wu, Proc. Natl. Acad. Sci. 105, 14262 (2008). 

14 F. Ma, W. Ji, J. Hu, Z.-Y. Lu, and T. Xiang, Phys. Rev. Lett. 102, 
177003 (2009). 



15 C. Fang, B. A. Bernevig, and J. Hu, EPL 86, 67005 (2009). 

16 O. J. Lipscombe, G. F. Chen, C. Fang, T. G. Pening, D. L. Aber- 
nathy, A. D. Christianson, T. Egami, N. Wang, J. Hu, and P. Dai, 
Phys. Rev. Lett. 106, 057004 (2011). 

17 J. Guo, S. Jin, G. Wang, S. Wang, K. Zhu, T. Zhou, M. He, and X. 
Chen, Phys. Rev. B 82, 180520 (2010). 

18 M.-H. Fang, H.-D. Wang, C.-H. Dong, Z.-J. Li, C.-M. Feng, J. 
Chen, and H. Q. Yuan, EuroPhys. Lett. 94, 27009 (201 1). 

19 A. Krzton-Maziopa, Z. Shermadini, E. Pomjakushina, V. Pom- 
jakushin, M. Bendele, A. Amato, R. Khasanov, H. Luetkens, and 
K. Conder, J. Phys.: Condens. Matter 23, 052203 (2011). 

20 Y. Mizuguchi, H. Takeyal, Y. Kawasaki, T. Ozaki, S. Tsuda, T. 
Yamaguchi, and Y. Takano, Appl. Phys. Lett. 98, 042511 (2011). 

21 D. M. Wang, J. B. He, T.-L. Xia, and G. F. Chen, Phys. Rev. B 83, 
132502 (2011). 

22 Y. Zhang, L. X. Yang, M. Xu, Z. R. Ye, F. Chen, C. He, H. C. Xu, 
J. Jiang, B. P. Xie, J. J. Ying, X. F. Wang, X. H. Chen, J. P. Hu, M. 
Matsunami, S. Kimura, and D. L. Feng, Nature Materials 10, 273 
(2011). 

23 T. Qian, X.-P. Wang, W.-C. Jin, P. Zhang, P. Richard, G. Xu, X. 
Dai, Z. Fang, J.-G. Guo, X.-L. Chen, and H. Ding, Phys. Rev. Lett. 
106, 187001 (2011). 

24 D. Mou, S. Liu, X. Jia, J. He, Y. Peng, L. Zhao, L. Yu, G. Liu, S. 
He, X. Dong, J. Zhang, H. Wang, C. Dong, M. Fang, X. Wang, 
Q. Peng, Z. Wang, S. Zhang, F. Yang, Z. Xu, C. Chen, and X. J. 
Zhou, Phys. Rev. Lett. 106, 107001 (2011). 

25 I. R. Shein and A. L. Ivanovskii, Phys. Lett. A 375, 1028 (2011). 

26 X.-W. Yan, M. Gao, Z.-Y. Lu, and T. Xiang, Phys. Rev. B 84, 
054502 (2011). 

27 C. Cao and J. Dai, Chin. Phys. Lett. 28, 057402 (201 1). 

28 L. Zhang and D. J. Singh, Phys. Rev. B 79, 094528 (2009). 

29 W. Bao, Q. Huang, G. F. Chen, M. A. Green, D. M. Wang, J. B. 
He, X. Q. Wang, and Y. Qiu, Chin. Phys. Lett. 28, 086104 (201 1). 

30 W. Bao, G. N. Li, Q. Huang, G. F. Chen, J. B. He, M. A. Green, 



9 



Y. Qiu, D. M. Wang, and J. L. Luo, arXiv: 1102.3674. 

31 F. Ye, S. Chi, W. Bao, X. F. Wang, J. J. Ying, X. H. Chen, H. D. 
Wang, C. H. Dong, and M. Fang, arXiv: 1 102.2882. 

32 R. Yu, J.-X. Zhu, and Q. Si, Phys. Rev. Lett. 106, 186401 (2011). 

33 Y. Zhou, D.-H. Xu, F.-C. Zhang, and W.-Q. Chen, Europhys. Lett. 
95, 17003 (2011). 

34 C. Cao and J. Dai, Phys. Rev. Lett. 107, 056401 (2011). 

35 Z. Wang, Y. J. Song, H. L. Shi, Z. W. Wang, Z. Chen, H. F. Tian, 
G. F. Chen, J. G. Guo, H. X. Yang, and J. Q. Li, Phys. Rev. B 83, 
140505 (2011). 



X.-W. Yan, M. Gao, Z.-Y. Lu, and T. Xiang, Phys. Rev B 83, 
233205 (2011). 

M. Wang, C. Fang, D.-X. Yao, G. Tan, L. W. Harriger, Y. Song, T. 
Netherton, C. Zhang, M. Wang, M. B. Stone, W. Tian, J. Hu, and 
P. Dai, arXiv: 1105.4675. 

Y.-Z. You, H. Yao and D.-H. Lee, Phys. Rev. B 84, 020406 (201 1). 
C. Fang, B. Xu, P. Dai, T. Xiang, and J. Hu, arXiv: 1103.4599. 
F. Lu and X. Dai, arXiv:l 103.5521. 



